Mon. Not. R. Astron. Soc. 416, 2068-2077 (201 1) Printed 28 January 2013 (MN IflgC style file v2.2) 



On the nature of the absorber in IRAS 09104+4109: the X-ray and 
mid-infrared view 
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ABSTRACT 

We present a long (sa 76 ks) Chandra observation of IRAS 09104+4109, a hyper-luminous 
galaxy, optically classified as a Type 2 AGN hosted in a cD galaxy in a cluster at z — 0.442. 
We also report on the results obtained by fitting its broad-band spectral energy distribution. 
The Compton-thick nature of this source (which has been often referred to as an 'archetype" 
of Compton-thick Type 2 quasars) was formerly claimed on the basis of its marginal detection 
in the PDS instrument onboard BeppoSAX, being then disputed using XMM-Newton data. 
Both Chandra analysis and optical/mid-IR spectral fitting are consistent with the presence of 
heavy («l-5 x 10 23 cm -2 ), but not extreme (Compton-thick) obscuration. However, using 
the mid-IR and the [O III] emission as proxies of the nuclear hard X-ray luminosity suggests 
the presence of heavier obscuration. The 54-month Swift BAT map shows excess hard X-ray 
emission likely related to a nearby (z=0.009) Type 2 AGN, close enough to IRAS 09104+4109 
to significantly enhance and contaminate its emission in the early BeppoSAX PDS data. 
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1 INTRODUCTION 

IRAS 09104+4109 (hereafter referred to as 109104) is a hyper- 
luminous infrared cD galaxy at z=0.442 residing in the core of a 
rich cluster of galaxies (Kleinmann et al. 1988). Its optical spec- 
trum is characterized by the presence of narrow emission lines, al- 
though broad Balmer and Mgll lines were observed in polarized 
light (Hines & Wills 1993; Tran, Cohen & Villar-Martin 2000). 
These remarkable results, based on optical observations, suggested 
the presence of a powerful quasar hidden in 109104. Subsequent 
X-ray observations confirmed this "picture". While BeppoSAX data 
below 10 keV were "contaminated" by the intra-cluster emission, 
the weak signal detected by the PDS instrument in the range 
~ 15-60 keV (in a 34 ks exposure) was interpreted by Frances- 
chini and collaborators (2000, hereafter F00) as the primary emis- 
sion of the buried quasar, emerging at very high energies because 
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of Compton-thick absorption (i.e., above ~ 10 24 cm -2 ; in the case 
of 109104, the estimated lower limit on the column density was 
5 x 10 24 cm -2 ). The detection of a strong iron Ka emission line, 
with equivalent width (EW) of ~l-2 keV - consistent with a re- 
flection scenario for the quasar emission in the 2-10 keV energy 
range - further supported the Compton-thick nature of 109104 in 
F00, and similar conclusions, though reached with limited counting 
statistics (but much higher spatial resolution), were drawn from a 
9.1 ks Chandra observation (Iwasawa, Fabian & Ettori 2001, here- 
after 101). 

Doubts on these conclusions were recently cast by Picon- 
celli et al. (2007, hereafter P07) on the basis of the analysis 
of a 14 ks XMM-Newton observation. In particular, although 
both a transmission- and a reflection-dominated model provided 
a good fit to the quasar emission, the iron line EW was lower 
(~ 400 eV) than expected in the case of a Compton-thick scenario 
(EW~ 1—2 keV), but it was higher than model predictions for 
a fewxlO 23 cm -2 obscuration (Ghisellini, Haardt & Matt 1994). 
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Furthermore, the source bolometric luminosity was well recovered 
using the 2-10 keV flux obtained in a transmission model (after re- 
moval of the significant contribution from the cluster, see P07) and 
assuming a reasonable X-ray bolometric correction (§4 in P07). 
According to these authors, the weak detection in the PDS instru- 
ment onboard BeppoSAX was explained by the presence of a heav- 
ily obscured, possibly Compton-thick galaxy (NGC 2782; Zhang et 
al. 2006) in the large (1.3° FWHM) field of view of the PDS. 

To date, the number of "certified" (i.e., secured by good- 
quality X-ray spectra and, possibly, high-energy detection) 
Compton-thick Active Galactic Nuclei (AGN) is still limited, be- 
ing settled around ~ 40 — 50 and mostly at low ( < 0.05) redshifts 
(see Comastri 2004 for a review on this topic, and Delia Ceca et al. 
2008 for un updated list of Compton-thick AGN). Large numbers 
of Compton-thick AGN candidates at higher redshifts have been re- 
cently claimed, but their identification as heavily obscured AGN is 
often derived from other, indirect arguments [e.g., optical emission- 
line strength, combined mid-infrared (mid-IR) and optical criteria, 
stacked X-ray emission; e.g., Alexander et al. 2008, 2011; Fiore 
et al. 2008, 2009; Lanzuisi et al. 2009; Bauer et al. 2010; Gilli et 
al. 2010, and references therein; Vignali et al. 2010] than direct, 
good-quality X-ray spectroscopy. 

Clearly, the lack of broad-band X-ray instruments within the 
same observatory (as in BeppoSAX) is currently limiting the detec- 
tion and spectral characterization of Compton-thick AGN; in this 
regard, the Integral-IBIS (Beckmann et al. 2009) and Swift BAT 
(Tueller et al. 2008) surveys have brought to the detection of addi- 
tional very local Compton-thick AGN candidates, but their confir- 
mation has always required follow-up observations with Chandra, 
XMM-Newton or Suzaku (e.g., Ueda et al. 2007; Eguchi et al. 2009; 
Comastri et al. 2010; Severgnini et al. 201 1 ; Burlon et al. 201 1). To 
overcome this observational problem will necessarily rely on the 
next-generation hard X-ray imager that will fly on NuSTAR. 

In this context, it is clear the role of our present investigation to 
assess the "true" nature of 109104, i.e., whether it is Compton thick 
or Compton thin. If Compton thick, it would be one of the few AGN 
of this class beyond the local Universe supported by moderately 
good-quality X-ray spectroscopy (for other moderate/high-redshift 
Compton-thick AGN, see Norman et al. 2002; Iwasawa et al. 2005; 
Comastri et al. 2011; Feruglio et al. 2011; Gilli et al. 2011). We 
tackle this issue through a multi-wavelength approach consisting 
of (a) the analysis of a long (~ 76 ks) and recent (2009) Chandra 
observation of 109104 and (b) the reconstruction of its optical and 
mid- and far-infrared (far-IR) spectral energy distribution (SED) in 
terms of stellar plus reprocessed AGN emission, using the code de- 
veloped by Fritz, Franceschini & Hatziminaoglou (2006, hereafter 
F06). While (a) will provide further clues on the best-fitting model 
characterizing the X-ray emission of 109104, (b) will allow us to 
estimate the geometry and coverage of the obscuring matter as well 
as its column density, which will be compared to that derived from 
X-rays. 

Finally, we use known correlations found in AGN to estimate 
the expected intrinsic 2-10 keV emission. In particular, we use the 
5.8|J.m, 12.3|J.m and [O III]5007A luminosities as a proxy of the nu- 
clear luminosity following the correlations reported by Lanzuisi et 
al. (2009) and Gandhi et al. (2009) for the mid-IR, and Mulchaey 
et al. (1994) and Heckman et al. (2005) for the optical. Further- 
more, we also compute the 2-10 keV over [Ne V]3426A flux ratio 
to possibly place further constraints on the amount of absorption. 

A A-cosmology with Hq=10 km s Mpc~' , £2^=0-3 and 
£2 A =0.7 (Spergel et al. 2003) is adopted. 



2 CHANDRA DATA 
2.1 X-ray data analysis 

109104 was observed on Jan 6 th 2009 for 76.16 ks with the Ad- 
vanced CCD Imaging Spectrometer (ACIS; Garmire et al. 2003) 
and the 13 CCD at the aimpoint. The observation (OBS JD=10445) 
was carried out in very faint mode. Since the source is in a mas- 
sive cooling-flow cluster (Fabian & Crawford 1995), showing a sig- 
nificant radial temperature gradient, the background subtraction is 
highly critical, despite the brightness of 109104. Therefore, to ex- 
tract spectra, we decided to adopt a strategy similar to that applied 
by 101 to the previous 9.1 ks Chandra observation, dating back to 
1999. We extracted the source spectrum from a 1 "-radius region 
(large enough to include ~ 85-90 per cent of the on-axis Chandra 
ACIS PSF) and the background from the surrounding l"-2" annu- 
lus, to possibly account for a significant fraction of the cool-core 
cluster thermal emission. The source extraction region was chosen 
from the X-ray image above 2 keV, where the contrast between 
109104 and the extended emission is maximized. At softer ener- 
gies, the source does not appear to be point-like, and its location 
seems consistent with the bright core of the cluster emission. Be- 
cause of its superior PSF, using Chandra data allows us to better 
evaluate the contribution due to the thermal emission - whose de- 
tailed analysis is beyond the purposes of this work - than in the 
XMM-Newton observation (P07). We note that after background 
removal, ~ 700 source counts are available for spectral fitting. The 
resulting spectrum has been grouped with a minimum of 20 counts 
per bin to apply the % 2 statistic and fitted with XSPEC V 12. 5.0 (Ar- 
naud 1996). 

The results obtained from this source spectral extraction has 
been checked - as reported in the following - using a different strat- 
egy, consisting of fitting the background contained in the l"-2" 
annulus, freezing its spectral parameters to their best-fitting val- 
ues, and then applying this background model (properly scaled to 
account for the factor 3 difference in area between source and back- 
ground) to the source spectral components (see 32.2) . 

Hereafter, errors are quoted at the 90 per cent confidence level 
for one interesting parameter (i.e., A% 2 = 2.71; Avni 1976); all 
spectral fits include absoiption due to the line-of-sight Galactic col- 
umn density of N H = 1.42 x 10 20 cirT 2 (Kalberla et al. 2005), and 
Anders & Grevesse (1989) abundances are assumed. 



2.2 X-ray spectral fitting 

To reproduce the Chandra data over the 0.5-7 keV energy range, 
we have adopted two models, similarly to the fitting strategy ap- 
plied by P07 to XMM-Newton data. The first model (hereafter re- 
ferred to as "transmission" model) comprises an absorbed power- 
law component plus an iron emission line. The second model (re- 
ferred to as "reflection" model) includes a reflection component 
(PEXRAV in XSPEC) and an iron emission feature. These models are 
reported as A and B in Table[T] where A2 differs from A1 because 
the powerlaw photon index is fixed to 1.8, as expected and typically 
found in case of AGN emission (e.g., Piconcelli et al. 2005, and 
references therein). Furthermore, both models include a soft com- 
ponent due to the residual (i.e., not totally accounted for by back- 
ground subtraction) cluster thermal emission. It is also possible that 
there might be a contribution to soft X-rays from photoionized gas 
at the inner nucleus, as suggested by 101; this hypothesis would be 
consistent with the hard-to-soft X-ray flux ratio (~ 4 per cent) ob- 
served in 109104, which is comparable to that typically measured 
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Table 1. Best-fitting spectral parameters for the Chandra spectrum of 109104. 



Model 

(1) 


kT 

(2) 


r 

(3) 


Nh Eica 
(4) (5) 


EW Ka 
(6) 


X 2 /dof 
(V) 






Source: R= 


=1"; Back: annulus 1"-2" 






A1 

A2 
B 


74 +5.73 

^•"'-1.58 

1 44+1.00 
1 -^-0.52 


70 +0 ' 71 
u./u_ 040 

1.8 (fr) 

1.8 (fr) 


2.59t\H 6.41±0.06 
4.29^82 6.41±0.06 
6.42±0.06 


4 25 +164 
4Z;, -I40 

316+ 210 
530+182 


30.0/39 
35.1/40 
53.2/40 



Source: R=1"; Back: annulus 1"-2" - Separate source and back fitting 



A1 3.81±t:lo OJ4 --°026 2.81+°^ 6.42±0.06 337+| 4 ? 26.4/39 
A2 4.02+ 4 j 8 1.8 (fr) 4.49^| 6.41+55;^ 274±118 30.8/40 
B 1.44^|j 1.8 (fr) 6.43±0.06 52.0/40 

Notes — The source counts are extracted from a circular region centered on 109104 of ra- 
dius 1", and the background from the surrounding annulus of size 1 "-2". The top part of the 
table refers to the "standard" case where the spectral parameters are obtained by subtracting 
the background to the source+background spectrum within XSPEC. The bottom part of the 
table reports the same spectral parameters when a different approach (fit to the background 
counts and subsequent fit to the source data, once the background components, properly 
normalized by the source area/background area, are frozen to their best-fitting values) is 
adopted. (1) X-ray modeling: Al: absorbed power-law model + iron line (transmission 
model); A2: the same, with T frozen to 1.8, as expected for AGN emission; B: reflection 
+ iron line (reflection model), with T frozen to 1.8. In both models, a thermal component 
has been included; (2) temperature of the gas in the cluster (fitted with a MEKAL model 
within XSPEC; (3) power-law photon index; (4) column density (in units of 10 23 cm -2 ); (5) 
rest-frame energy of the iron Ka line (in keV), whose width is fixed to 10 eV; (6) rest-frame 
equivalent width of the iron line (in eV); (7) x 2 /de grees of freedom (dof). 



in Sey 2 galaxies (e.g., Matt et al. 2004; Piconcelli et al. 2008, and 
references therein). However, a detailed analysis of the soft X-ray 
emission, limited by the low number of counts below 2 keV (~ 200) 
and the lack of deep X-ray grating data, is beyond the purposes of 
this paper, since the main goal of our work consists of a proper 
characterization of the hard (>2 keV) emission. 

In the upper part of Table [Tj we report the X-ray spectral 
results for both models using a "standard" background subtrac- 
tion within XSPEC, where the source+background spectrum is ex- 
tracted from a circular region of radius 1" centered on the posi- 
tion of 109104, and the background spectrum from a surround- 
ing annulus, as anticipated in 92.11 The lower part of Table Q] re- 
ports the same models using a different background subtraction 
strategy; at first, the background, using its own response matri- 
ces, is fitted with a thermal component plus a low-normalization 
powerlaw component. The best-fitting parameters are then frozen, 
and the source+background spectrum is fitted using these compo- 
nents (multiplied by a normalization constant fixed to 0.33, i.e., the 
source/background area ratio) and those ascribed to the nucleus of 
109104. As evident from Table [T] and already reported above, there 
is still need for a thermal component, possibly due to the residuals 
at soft X-ray energies obtained from the best-fitting of the back- 
ground spectrum. By using the latter background subtraction pro- 
cedure, errors on the parameters are typically lower. 

Whatever the method adopted to account for background- 
subtraction is, the transmission model is preferred in terms of qual- 
ity of the fit (see column (7) in TableQJ. In this model, the column 
density is Njj~ (1 — 5) x 10 23 cm~ 2 , closely resembling the ab- 
sorption values found by P07. The iron line, compatible with being 
neutral or mildly ionized, has an EW~ 300 — 400 eV (in the rest 
frame of the source), is consistent, within errors, with the value 
found by P07 (and similarly for the line flux) and with that ob- 



tained by 101 using lower exposure Chandra data. We note that 
these large EW values, though not extreme, are slightly higher than 
expected in the case of fluorescence iron line produced by the ob- 
scuring matter (Ghisellini et al. 1994^] but appear consistent with 
previous XMM-Newton and Suzaku observations of Compton-thin 
Seyfert 2 galaxies (e.g., Guainazzi, Matt & Perola 2005; Fukazawa 
et al. 2011). Given the large stellar mass estimated via SED fit- 
ting (which is described in 5j4), ~ 4.8 x 10 11 M Q , and the 
known mass-metallicity relation (Lequeux et al. 1979; Savaglio et 
al. 2005), we may expect some metal (more specifically, iron) en- 
hancement in 109104, hence an increased strength of the Fe emis- 
sion line, although it is not straightforward to quantify such effect. 
For instance, NGC 4388, a local Sey 2 galaxy with a column den- 
sity comparable to that of 109104, shows a similarly large iron Ka 
EW (Iwasawa et al. 2003; see also Beckmann et al. 2004 and Shirai 
et al. 2008 for significant EW differences due to continuum vari- 
ability), although it is not hosted by a massive galaxy. On the other 
hand, we note that Mrk 231 shows a relatively low EW for the 
iron Ka line, despite being a Compton-thick AGN (Braito et al. 
2004). More convincing hints for a Compton-thin absorption to- 
wards 109104 are likely provided by the high-energy constraints 
(see HJl. 

The 90 per cent upper limit on the EW for He-like and H- 
like iron line at 6.7 keV and 6.97 keV is ~ 85 eV. Although the 
powerlaw is characterized by an apparently flat photon index, the 
solution with a more physically acceptable T = 1.8 still provides 
an adequate fit to the data (see TableQ}. The source spectrum, fitted 



The predicted EW for the iron line depends on the inclination angle of 
the torus, although this dependence is not strong given the column density 
range derived from X-ray spectral fitting of 109104 data. 
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Figure 1. Chandra datapoints of 109104 fitted with an absorbed powerlaw 
plus an iron line (transmission model; A1 model in the top part of Table[T] 
and thermal emission, ascribed to the cluster. 

using a transmission model (including a neutral iron line), is shown 
in Fig.Q] 

The 2-10 keV flu£|is « (3.8-4.2) x 10~ B erg cirr 2 s _1 , 
corresponding to an intrinsic (i.e., absorption-corrected in the case 
of models A1 /A2) AGN luminosity (placing the contribution due to 
the residual thermal emission to 0) of ~ (2.7 — 3.1) x 10 44 erg s 
(where the range of values is ascribed to the differences in the 
method adopted for background subtraction). 

Although the low-energy (below ~ 6 keV) range provided 
by Chandra is apparently not well suited for placing significant 
constraints on the reflection component, we modeled the data us- 
ing a Compton reflection component from neutral matter (PEXRAV 
model in XSPEC; Magdziarz & Zdziarski 1995); we adopted cos 
i=0.9 as inclination angle for the reflection component, which is 
equivalent to assuming a rather face-on geometry for the reflecting 
material, likely to be associated with the inner walls of the obscur- 
ing torus; we chose a cut-off energy of 200 keV for the primary 
continuum, and metal abundances were fixed to the solar value. We 
also fixed the photon index of the continuum to 1.8 to facilitate 
the spectral fitting with the reflection model over the limited en- 
ergy bandpass provided by Chandra and get a reasonable constraint 
on the strength of such component. The Chandra data require a 
reflection-dominated spectrum, where the primary powerlaw nor- 
malization is formally consistent with 0. From a statistical point of 
view, the reflection model provides a worse description of the X-ray 
data, regardless of the adopted method to subtract the background 
(see Table QJ. The iron Ka emission line EW is ~ 400 — 500 eV, 
which is not easily reconcilable with what expected in a typical 
reflection scenario. 

We note that the photon index derived by our spectral fittings 
in case of a transmission model (Al) is flat. This result can suggest 

2 Since the currently used version of PSEXTRACT, the tool adopted to ex- 
tract the X-ray spectrum from Chandra data, does not take into account 
the correction for the encircled energy fraction enclosed in the source ex- 
traction region, we compute average "correction factors" to be applied to 
the fluxes using the Chandra Proposer's Guide. At the on-axis position of 
our target, these factors, included in the flux and luminosity values reported 
hereafter, are Rt 10 and rj18 per cent for the 0.5-2 keV and 2-10 keV bands, 
respectively. 



that an additional reflection component may be present, causing the 
observed X-ray photon index to be flatter than expected from AGN 
emission (e.g., Piconcelli et al. 2005); however, a "mixed" absorp- 
tion and reflection model is not able to provide any improvement to 
the previous spectral solutions, calling for higher quality data over 
a possibly more extended energy interval. 

Overall, it seems that any residual cluster emission has been 
taken into account properly by our procedure of background re- 
moval and fitting, thanks to the good on-axis PSF provided by 
Chandra. For instance, the ionized iron line components related 
to extended (cluster) emission may be at the origin of the large iron 
line EW previously reported (F00). 

The Chandra data presented in this paper provide consistent 
results with those published by P07; although the photon statistics 
is higher in the XMM-Newton observation (which also guarantees 
data over a broader energy range), Chandra data allow us for a bet- 
ter treatment of the contamination from the cluster in virtue of its 
better PSF. We find that ~ 90 per cent of the flux is ascribed to 
the active nucleus over the 2-10 keV range, to be compared with 
~ 30-35 per cent derived by P07, where the source counts were ex- 
tracted from a much larger [37" (pn) - 40"-radius (MOS)] region. 



3 HIGH-ENERGY CONSTRAINTS FROM THE 
54-MONTH SWIFT BAT DATA 

We used the 54-month Swift BAT data, processed using the 
BATjMAGER software (Segreto et al. 2010), to constrain the high- 
energy emission of 109104. The source is not detected in any of 
the BAT sub-bands, the 3-fJ upper limits to the source X-ray emis- 
sion (derived using a Crab-like spectrum) being ~ 1.9 x 10 -12 , 
fa3.3x 10~ 12 ,and^4.8x 10~ 12 erg cm" 2 s" 1 in the 15-30 keV, 
15-70 keV, and 15-150 keV bands, respectively. The lowest energy 
upper limit provided by BAT is consistent, within uncertainties, 
with the flux derived in the 20-30 keV energy range by F00 using 
BeppoSAX data (sa 2.5 x 10~ 12 erg cm~ 2 s -1 ). The flux limit in 
the BAT map at the position of 109104 (but see below for concerns 
about possible contamination) is too loose to place any constrain 
on the two X-ray models presented in 32.21 

As shown in Fig. [2] there is "excess" emission (at the ~ 4a 
level, i.e., just below the detection threshold adopted for BAT maps) 
close (at ~ 17') to our source of interest, although such emission, 
at first glance, appears related to another object. It is reasonable to 
believe that the upper limit to the hard X-ray emission of 109 104 re- 
ported above is actually contaminated, at least in part, by the nearby 
source's PSF wings and, as such, must be intended as very con- 
servative and poorly indicative of the effective high-energy emis- 
sion of 109104. This finding provides support to the hypothesis of 
a Compton-thin absorber towards 109104. 

Using the NASA Extragalactic Database (NED) and literature, 
we have attempted to assign a known counterpart to the source 
closest to 109104. Taking into account the positional uncertainty 
in the BAT map (~ 6', corresponding to the size of the circles in 
Fig. [2}, at the position of the BAT hard X-ray emission we found 
NGC 2785, a disky galaxy at z = 0.009, whose Sloan Digital Sky 
Survey (SDSS) spectrum is characterized by a very red continuum 
and a large Hoc/HP ratio (~ 20.9 at face value, or ~ 12.5 taking 
into account the correction for the average host galaxy absorption 
component), being highly suggestive of a heavily extincted object. 
The "corrected" Hoc/Hp ratio is consistent with the value reported 
by Veilleux et al. (1995) using a different optical spectrum. Un- 
der reasonable assumptions, this line flux ratio translates into a 
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Figure 2. Swift BAT 54-month significance map (15-30 keV band) of 
the area around 109104, clearly showing hard X-ray emission close but 
likely not associated to the z = 0.442 quasar. The position of NGC 2782, 
previously claimed to contaminate the PDS flux of 109104 (P07), is 
also shown, along with the positions of NGC 2785 (black circle) and 
SDSS J091548.16+410254.23; see text for details. All circles have radii 
of 6', broadly corresponding to the positional uncertainty provided by BAT 
for weak sources. The color scale indicates the intensity level of the pixels 
in the map. 



E(B-V)»1.5-1.9. Furthermore, the Hoc/[NII] and Hp/[OIII] line 
ratios are consistent with this source being a Type 2 AGN. The as- 
sociation of the BAT excess with this source appears as the most 
likely possibility, as also suggested by their good spatial coinci- 
dence (see Fig.O. At present, the only X-ray constraint to the emis- 
sion of NGC 2785 below 10 keV (over a broader energy range than 
that allowed by the shallow ROSAT All Sky Survey, RASS, where 
the source is not detected) can be derived by an ASCA GIS ob- 
servation, where the source, located at a large off-axis angle from 
the nominal pointing, is not detected. Assuming a simple power- 
law model with T=1.8, the derived 3a upper limit to its 2-10 keV 
flux is ~ 1.2 x 10~ 13 erg cm~ 2 s _1 . Using the de-reddened [O III] 
emission to estimate the intrinsic 2-10 keV flux (adopting the cor- 
relation found by Mulchaey et al. 1994), we derive an X-ray flux of 
~2.9x 10~" erg cm~ 2 s . At face value, the difference between 
the upper limit to the hard X-ray flux derived from ASCA data and 
that estimated from the [O III] intensity can be explained assuming 
an absorber with a column density of at least ~ 2.6 x 10 24 cm -2 . 
However, the upper limits to the X-ray emission of NGC 2785 
above 10 keV derived from the BAT instrument are well below the 
extrapolation of the absorbed powerlaw, thus casting some doubts 
on either the X-ray flux derived from the [O III] measurement or the 
adopted model (or both). Pointed and sensitive X-ray observations 
of NGC 2785 are required to shed light on all these issues. 

At a close position, we found another possible (but less 
likely) optical counterpart, a broad-line SDSS quasar at z = 2.32, 
SDSS J091548. 16+410254.23. This quasar has a tentative detec- 
tion (^ 2.4o) in the RASS at w 19". However, from the BAT map, 
it seems clear that this source does not provide a significant contri- 
bution to the observed signal. 

Although the association of the high-energy emission with 
NGC 2785 requires further support, we note that BAT data have 
clearly shown that contamination in the early PDS data is a valu- 



able hypothesis. As a consequence, serious doubts are cast on the 
Compton-thick nature of 109104, since it was mostly motivated by 
the admittedly marginal (~ 2.5a) detection in the PDS instrument 
and the large equivalent width of the iron Ka line (F00), which 
most probably suffered from the contribution of ionized iron fea- 
tures associated to the thermal emission in BeppoSAX data. 

Furthermore, from the BAT map it seems clear that NGC 2782 
(z = 0.008), the Seyfert 2 galaxy thought to contaminate the hard 
X-ray emission of 109104 in the PDS (P07), has marginal emission 
in the hard X-rays (see Fig.|2j, and its distance from the z = 0.442 
quasar (~ 50') is beyond any possible contamination in the BAT 
map. 



4 SPECTRAL ENERGY DISTRIBUTION: 

MULTI- WAVELENGTH DATA AND FITTING 

To provide a better understanding of the matter responsible for the 
X-ray absorption towards 109104, we have investigated the source 
broad-band properties using the data available in literature. The sci- 
entific goal consists of using the reprocessed emission to constrain 
the physical/geometrical parameters of the absorber (the "postu- 
lated" torus in AGN unification schemes) and estimate the bolo- 
metric luminosity of the AGN due to accretion of matter onto the 
black hole. This study requires that all the possible source com- 
ponents emitting in the UV, optical, near-IR, mid-IR, and far-IR 
bands are properly taken into account in terms of stellar (galaxy) 
plus nuclear (AGN) components. The broad-band analysis has been 
carried out by using most of the photometric datapoints available 
in Table B.4 of Ruiz et al. (2010). In particular, in our SED fit- 
ting analysis, we used the XMM-Newton Optical Monitor (OM) 
data (filters UVW2 and UVW1), the SDSS magnitudes (excluding 
the /-band magnitude, given its contamination from strong [O III] 
emission), the published Spitzer IRAC data at 3.6(xm and 5.8|J.m, 
and MIPS data at 24ii.m. The source was also detected by IRAS at 
all wavelength^ up to 60ii.m, while at longer wavelengths (IRAS 
at 100ii.m, SCUBA at 450ii.m and 850(xm), only upper limits are 
available. We remind the reader of the non-simultaneity for most 
of the data used in our analysis. This issue, common to most in- 
vestigations dealing with SED fitting, likely plays a minor role, 
given the long timescales for variability expected in case of repro- 
cessed emission by dust. We also note that our photometric points 
include both nuclear (i.e., AGN-related) and galaxy emission, with 
both components being accounted for by the adopted broad-band 
SED-fitting procedure, as described in the following. Finally, as a- 
posteriori check, we have verified that the broad-band modeling 
adopted to reproduce the emission from 109104 is consistent, in the 
mid-IR, with the Spitzer IRS spectrum available in literature (dis- 
cussed in 94. 1 \ . 

Data have been modeled using the code developed by F06, 
which assumes that the dust, responsible for the reprocessing of the 
nuclear radiation, consists of graphite and silicate grains distributed 
in a "flared disc", "smooth" geometry. Recent high-resolution, in- 
terferometric mid-IR observations of nearby AGN (e.g., Jaffe et 
al. 2004; Tristram et al. 2007; Meisenheimer et al. 2007; Tris- 
tram & Schartmann 2011) have confirmed the presence of a ge- 
ometrically thick, torus-like dust distribution on pc-scales; how- 
ever, data strongly suggest that this torus is irregular or "clumpy". 

3 For the flux densities at the IRAS wavelengths, we refer to the association 
of 109104 with the IRAS Faint Source Catalog object F09105+4108. 
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Table 2. Best-fitting SED parameters using a "smooth" torus model. 
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Notes — For a full description of the parameters, see Fritz et al. (2006). 
(1) Angle between the line of sight and the equatorial plane of the torus (in 
degrees); (2) optical depth at 9.7 u.m; (3) full covering angle of the torus (in 
degrees); (4) E(B — V) of the stellar component; (5) "merit function" (sort 
of x 2 ), used to evaluate the goodness of the fit (see text for details). 

X-ray variability on short time-scales provides further support to 
such clumpy-absorption scenario (Risaliti et al. 2007, 2010). In this 
regard, many codes have been developed over the last decade to 
deal with clumpy dust distributions (e.g., Nenkova, Ivezic & Elitzur 
2002; Nenkova et al. 2008a,b, N08a,b hereafter; Schartmann et al. 
2008; Honig & Kishimoto 2010). We note, however, that fitting 
broad-band SED data, based on sparse photometric datapoints, with 
either smooth or clumpy model is not supposed to provide signif- 
icantly different results (see Dullemond & van Bemmel 2005) in 
terms of derived AGN bolometric luminosity. In the following, we 
will mainly focus on the best-fitting parameters obtained by adopt- 
ing the smooth dust distribution (F06); however, a discussion based 
on the broad-band data fitting using the clumpy code described by 
N08a,b will also be provided. 

The SED fitting is based on a multi-component analysis, 
which has been adopted often in recent literature for the analysis of 
obscured AGN (e.g., Vignali et al. 2009; Hatziminaoglou, Fritz & 
Jarrett 2009; Pozzi et al. 2010). The observed UV to far-IR SED of 
109104 has been de-composed in three distinct components: stars, 
having the bulk of the emission in the optical/near-IR; hot dust, 
mainly heated by UV/optical emission due to gas accreting onto 
the supermassive black hole and whose emission peaks somewhere 
between a few and a few tens of microns; cold dust, principally 
heated by star formation. 

The stellar component has been included using a set of Sim- 
ple Stellar Population (SSP) spectra of solar metallicity and ages 
ranging from ~ 1 Myr to ~7.3 Gyr, which corresponds to the time 
elapsed between z=4 (the redshift assumed for the stars to form) and 
z = 0.442 in the adopted cosmology. A common value of extinction 
is applied to stars of all ages, and a Calzetti et al. (2000) attenua- 
tion law has been adopted (Ry = 4.05). Assuming a higher redshift 
for the event of star formation does not provide significantly differ- 
ent results, while a Milky Way extinction curve produces a slightly 
worse fit than the one reported here (see Table |2j- 

The AGN emission, dominating in the mid-IR range, com- 
prises a thermal component due to reprocessing of the nuclear AGN 
emission, a scattering component and a disc component. In partic- 
ular, we used an extended version of the model grid described by 
F06. 

Given the presence of data above 24(jm, also an emission 
component coming from colder, diffuse dust, likely heated by star- 
formation processes, has been included in the fitting procedure, and 
it is represented by templates of known starburst galaxies (see Vig- 
nali et al. 2009 for details). 

Similarly to Vignali et al. (2009) - to which we refer for fur- 
ther details on the adopted procedure and its limitations - we eval- 
uate the goodness of the fit using a "merit function", defined as 

MF = £^?j f^^) /N , where N is the number of observed 
datapoints, M,- and O; are the model and observed flux densities, 



T 



1 r 




Observed — frame Wavelength (/^m) 



Figure 3. Observed- frame SED of 109104, from the XMM-Newton Optical 
Monitor (OM) ultraviolet datapoints to the SCUBA 850|J.m datapoint. 3-a 
upper limits refer to non-detections (IRAS 100u.m and both SCUBA wave- 
lengths). The total SED (thick black line) is the summed contribution of 
the stellar (green dotted line) and AGN components (blue shaded line). The 
continuous red line over the wavelength range ss 5 — 38u.m represents the 
IRS spectrum, without any normalization applied. We note that the best- 
fitting solution obtained using the F06 model is consistent with the IRS 
spectrum in terms of normalization and reproduction of the 9.7u.m silicate 
absorption feature. 

respectively, of the i — th photometric band, and o,- is the corre- 
sponding observed error. 

From the SED fitting, there is no need for invoking cold dust 
ascribed to starburst emission, as suggested by the absence of de- 
tected emission in the far-IR above 60|i.m (Fig.[3j- The upper limits 
at longer wavelengths appear fully consistent with AGN emission 
(see also Rowan-Robinson 2000 for similar conclusions). In the 
UV/optical band, the stellar component has a moderate extinction 
[E(B-V)=0.16]. 

The best-fitting solution (Fig. [3} in terms of MF (see Table[2} 
provides a torus full covering angle of 140°, corresponding to a 
covering factor of ~ 90 per cent. The half opening angle, 20° 
((180-140/2)), is larger than the one suggested by the presence of 
a well collimated radio jet (Kleinmann et al. 1988; Hines & Wills 
1993) and a narrow optical ionization cone (Crawford & Vanderri- 
est 1996); see Fig. 3 of Taniguchi et al. (1997). We note, however, 
that none of the torus models in the grid used in this paper has 
a covering angle larger than 140°; to limit the degeneracies in the 
model parameters, only covering angles of 60, 100, and 140 degrees 
are actually used (see Fritz et al. 2006). These values refer to three 
possible "classes" of geometric solutions for the torus: a disk-like 
torus, a moderate-coverage torus, and one covering most of the nu- 
clear source. Our best-fitting SED delineates a configuration where 
a significant fraction of the innermost regions of 109104 is hidden 
by the torus; in these terms, its covering angle can be considered 
consistent with indications obtained with previous observations. 

The extinction at 9.7 (jm along the equatorial plan^3 is rela- 
tively high (x(9.7)=3) but not extreme, as one might expect in case 



Since the angle between the line of sight and the equatorial plane is zero, 
the extinction along the line of sight and that along the equatorial plane co- 
incides. Angles different from would imply an extinction in the direction 
of the observer lower than the equatorial one. 
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of Compton-thick absorption, and corresponds to a column density 
of ~ 1.1 x 10 23 cm~ 2 (assuming a Galactic gas-to-dust ratio con- 
version). This can be considered broadly consistent with those de- 
rived from X-ray spectral analyses [from both XMM-Newton (P07) 
and Chandra data; see 92.21 . Given the discrete values for the torus 
model grid parameters, it is not straightforward to associate un- 
certainties to the reported values but it is possible to evaluate how 
much solid the best-fitting solution is. Analysis of the degeneracies 
in smooth torus parameters (see Pozzi et al. 2010 for details) and 
corresponding uncertainties indicates that the value of the cover- 
ing angle of the torus (140°) is well determined, and x(9.7) along 
the equatorial plane reaches at most the value of 6 (90 per cent 
confidence level), which means a doubling of the column den- 
sity of the absorber. In all the solutions provided by F06 code for 
109104, an absorber with a corresponding column density of at least 
~ 10 23 cm -2 (but never extreme) is required by the fit. 

Overall, the best-fitting solution for the stellar plus nuclear 
emission of 109104 (Table O obtained using, for the reprocessed 
nuclear component, the smooth dust distribution (F06), is able to 
reproduce most of the observed photometric datapoints as well as 
the shape of the IRS spectrum and the 9.7 \im silicate absorption 
feature, as discussed in 94.11 The resulting accretion-related bolo- 
metric luminosity is ~ (2.3 — 3.8) x 10 47 erg s (where the range 
provides approximately the 90 per cent uncertainty on the luminos- 
ity). 

Given the relatively large number of photometric data avail- 
able and the presence of mid-IR spectroscopic coverage for 109104, 
this source can be considered a good target for a comparison of 
smooth vs. clumpy dust distribution; however, for an exhaustive 
comparison on a large AGN sample, we refer to Feltre et al. (201 1). 
Similarly to the SED fitting described above, we used a stellar plus 
AGN emission to account for the broad-band emission of 109104. 
In particular, for the latter component we used the clumpy dust 
models from N08a,b; according to these models, the dust is dis- 
tributed in clouds, hence there is always a finite probability for 
an unobscured (i.e., direct) view of the AGN, irrespective of the 
viewing angle (e.g., Nikutta, Elitzur & Lacy 2009). Among the 
latest publicly available grid models (end of September 2010), 
we adopted those which are considered best suited for Type 2 
AGN, where the direct component is likely minor or negligible (R. 
Nikutta, private communication). This choice is also motivated by 
the necessity of considering the contribution, at optical and near-IR 
wavelengths, of the host galaxy stellar component, which is visible 
in optical images. Including both stellar and direct nuclear com- 
ponents would have provided a over-estimation of the emission at 
wavelengths short-wards a few microns by the model. 

The clumpy model provides a good fit to the observed data 
in terms of data- to-model residuals (MF~4.4). However, we note 
two major problems in dealing with this solution (hence, with the 
parameters derived for the clumpy torus): (a) the model is not able 
to reproduce the near-IR (3.6|J.m) data, irrespective of the assumed 
stellar population fitting the optical data. As extensively discussed 
by Mor, Netzer & Elitzur (2009) and, more recently, by Deo et al. 
(201 1), this emission, not properly accounted for by Nenkova et al. 
grid models, needs an extra hot-dust component (parameterized by 
a blackbody with temperature of ~ 1400 K), possibly associated 
with graphite grains (which, given a higher sublimation temper- 
ature than silicate grains, can survive closer to the black hole; see 
Mor & Trakhtenbrot 201 1). At face value, this result may cast some 
doubts on the parameters derived in the mid-IR in case of clumpy 
models. In fact, the extra model component appears hard to justify, 
since it should be already accounted for by the radiative transfer 



treatment, as in smooth torus models. We note, however, that using 
the Draine & Lee (1984) optical properties for the silicate compo- 
nent (as in Nenkova et al. 2002) instead of the Ossenkopf, Henning 
& Mathis (OHM; 1992), which have been assumed in N08a,b and 
recent works (Nikutta et al. 2009), the discrepancy between model 
and data in the near-IR decreases (see Fig. |4). The second major 
concern is that (b) the best-fitting clumpy solution to the photomet- 
ric data is not able to reproduce the shallow 9.7 \im absorption fea- 
ture due to silicates in IRS spectrum. This result, which will be dis- 
cussed further in 94.11 is likely due to the fact that such SED-fitting 
solution requires the presence of many clouds along the equatorial 
plane. Their combined contribution manifests itself in a deep sili- 
cate feature, which is actually not observed as such in IRS spectrum 
(see Fig.[4]l. 

Since smooth-torus models are able to account, in a self- 
consistent way, all of the major components related to direct, scat- 
tered and reprocessed AGN emission, we prefer this solution for 
109104. However, we caution the reader against over-interpreting 
our simplistic comparison between the smooth and the clumpy so- 
lutions. Extended and more exhaustive checks need to be carried 
out on a sizable sample of AGN (Feltre et al. 201 1) before drawing 
general conclusions on the capability, by either model, to account 
for the reprocessed nuclear emission. 

For what concerns the smooth torus solution, we note that the 
SED fitting presented by Taniguchi et al. (1997) suggested the ad- 
ditional presence of cold dust distributed in a much larger outer 
torus, not included in the radiative transfer modeling. This dust 
layer would increase the silicate absorption at 9.7 micron, which 
has been found to be relatively shallow in IRS data. About the other 
"inner" torus parameters, their best-fitting solution implies a large 
covering factor (0.98, not too different from our solution; see dis- 
cussion in §4$ and a larger inclination with respect to the equatorial 
plane (30°) than ours. 

Finally, we may compare the outcome of our SED fitting in 
terms of gas mass with results from literature. Assuming a typical 
dust-to-gas mass ratio of 0.01, the dust mass obtained by our SED 
fitting implies a large gas mass (~ 5 x 10 10 Mq), a factor ~ 2 — 3 
higher than the upper limit derived by Evans et al. (1998) through 
CO measurement. However, both estimates have large uncertain- 
ties, mainly related to the underlying assumptions (e.g., significant 
variations in the dust-to-gas mass ratio across different galaxy pop- 
ulations have been recently reported by Santini et al. 2010). 

4.1 IRS spectrum of 109104 

Our own analysis of the Spitzer IRS spectrum, covering a wave- 
length range ~ 5 — 38iim in the observed frame, confirms the re- 
sults reported by Sargsyan et al. (2008; see also Deo et al. 2009); 
besides the 9.7iim silicate feature in absorption, suggestive of ob- 
scuration, the IRS spectrum is featureless, except for the very shal- 
low silicate absorption feature, and entirely consistent with the ob- 
served mid-IR datapoints (see Fig.[3]l. The constraints available for 
the typically most prominent PAH features, not observed in 109104 
(see also Sargsyan et al. 2008), are consistent with the lack of a 
significant starburst component, as found from both the SED fitting 
analysis (see §4$ and using IRAS (Rowan-Robinson 2000) and ISO 
data (Taniguchi et al. 1997). For completeness, we note that in case 
of a strong radiation field (as in the case of our target), the PAHs 
can hardly survive, if not properly screened or produced in regions 
far away from the active nucleus. 

We may use the IRS spectrum to make some further consider- 
ations about the smooth vs. clumpy solutions for the matter respon- 
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Figure 4. Observed-frame SED of 109104 obtained by fitting the same 
photometric points as in Fig. [3] The stellar emission (green dotted line) is 
accounted for using a set of SSPs, while the AGN emission is reproduced 
using a clumpy dust distribution (blue shaded line). The total SED is marked 
by the thick black line, while the IRS spectrum is plotted as a continuous 
red line. The Draine & Lee (1984) optical properties for the silicates have 
been assumed here (see text for details). 

sible for the reprocessing of the nuclear emission in the mid-IR 
described above. For what concerns the smooth-torus solution, we 
note that fitting the IRS spectrum alone (without taking into ac- 
count any stellar component, because of the IRS wavelength cover- 
age) provides, as best-fitting solution, exactly the same found to re- 
produce the source photometry (see Table|2}- Both the shape of the 
continuum and the 9.7 (im silicate absorption feature are accounted 
for reasonably well (see Fig. |3}. About the clumpy solution, using 
the Draine & Lee (1984) optical properties for the silicates pro- 
duces a good fit to the IRS continuum and 9.7(j.m silicate feature 
which, on the other side, is never properly reproduced by clumpy 
models with OHM silicate dust composition. This is opposite to 
what has recently been reported in literature (e.g., Nikutta et al. 
2009). Overall, our results on clumpy tori applied to 109104 in- 
dicate that the OHM composition may be rejected, being unable to 
reproduce the near-IR datapoints and the silicate absoiption feature. 
Clearly, this result, based on a single source, needs to be tested on 
larger samples. We finally note that in neither case the best-fitting 
clumpy solution to the IRS spectrum is the same one obtained by 
fitting the photometric data. 



5 DISCUSSION 

Two of the evidences that brought F00 towards an interpretation of 
109104 as a Compton-thick AGN, i.e., the hard X-ray emission and 
the large iron Koc line EW, have been discussed both by P07 us- 
ing the XMM-A'ewton data and in this paper using Chandra data. 
While the emission above 10 keV detected by the PDS instrument 
onboard BeppoSAX is very likely associated to a local, obscured 
AGN (see ^3}, the large iron emission line EW observed by Bep- 
poSAX was probably originated by a blend of emission features due 
to the AGN and the cluster. The latter emission has been properly 
accounted for by the analysis adopted in this paper thanks to the 
higher Chandra spatial resolution with respect to previous X-ray 
data which has limited the contamination by extended emission 



(see j )2,2t . The analysis of Chandra data, coupled with the SED re- 
sults (under the assumption of a smooth dust distribution), suggests 
that Compton-thin absorption is a likely solution for the absorber 
close to 109104. However, the upper energy boundary covered with 
sufficient statistics by Chandra is ~ 6 keV (observed frame); an 
extended energy range, possibly coupled with simultaneous data 
above 10 keV, would be useful to further support this transmission 
scenario. To provide a broader view of the possible solutions for 
such absorber, in the following we use known relations to estimate 
the intrinsic 2-10 keV luminosity of the source (to be compared 
with the observed luminosity), once the cluster emission has been 
taken into account. The following relations are characterized by 
significant dispersions; using them for a single object may suffer 
from significant limitations, therefore any conclusion drawn in the 
following should be treated with care. 

Given the AGN dominance in the mid-IR, the IRS spectrum 
of 109104 is used to derive the rest-frame 12.3 pm flux density and, 
adopting the Gandhi et al. (2009) correlation (which, obtained from 
a sample of 42 Seyfert galaxies with near diffraction-limited mid- 
IR imaging, links the mid-IR luminosity to the 2-10 keV luminos- 
ity; see, for a recent application, Vignali et al. 2010), we are able to 
estimate the expected X-ray emission for this source. On the basis 
of the SED fitting results (see 5j4), we do not expect a significant 
contribution from star formation at these wavelengths. 

The predicted AGN 2-10 keV luminosity 
(~ 8.9 x 10 45 erg s _1 ) is higher than the measured value 
(i.e., not corrected for the absorption) in case of a transmission 
scenario (~ 1.2— 1.3 x 10 44 ergs -1 ). The predicted/observed 
luminosity ratio is large and suggestive of Compton-thick obscura- 
tion. However, if we use the 5.8pm luminosity and the correlation 
reported in Lanzuisi et al. (2009), which has been obtained for 
a sample of quasars (i.e., at luminosities comparable to that of 
109104), albeit of limited size, we find that the expected 2-10 keV 
luminosity is ~ 10 45 erg s , much closer to the observed source 
luminosity but still suggesting heavy obscuration. 

The hard X-ray emission of 109104 can be estimated also us- 
ing [O IIl]-based correlations, under the assumption, motivated by 
unified schemes for AGN, that the [O III]5007A line is a good 
proxy of the nuclear emission. Using the value for the [O III] lu- 
minosity of 109104 taken from Reyes et al. (2008), not corrected 
for extinction within the narrow-line region (NLR), we obtain a 
2-10 keV luminosity of 5.2xl0 45 erg s" 1 and 3.5xl0 45 erg s" 1 
using the correlation, valid for Type 1 AGN, of Mulchaey et al. 
(1994) and Heckman et al. (2005), respectively. These values lie 
in between those obtained from mid-IR estimators and discussed 
above. 

As extensively described in Vignali et al. (2010), these cor- 
relations may suffer from the still limited number of objects used 
to derive them and, in particular, to the small number of AGN at 
high X-ray luminosity, i.e., in the quasar regime (indicatively above 
10 44 erg s _1 , where also 109104 lies). For these reasons, once more 
we advise against over-inteipreting our estimates of the intrinsic 
nuclear luminosity via the methods reported here, given the overall 
large uncertainties which characterize these correlations and their 
application to a single object. However, these results may indicate 
that some extra absorption with respect to that currently measured 
by Chandra might be present in 109104. 

Finally, we can measure the ratio between the AGN rest-frame 
2-10 keV flux prior to absorption correction and the [Ne V] emis- 
sion line (as available from the SDSS spectrum), and compare the 
resulting value (~ 80-90) to the recent findings published by Gilli 
et al. (2010). This ratio is expected to decrease for increasing col- 
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umn density, since the X-ray emission suffers from photo-electric 
absorption, while the [Ne V] emission, produced in the NLR, can be 
considered a proxy of the intrinsic nuclear power, although, given 
its rest-frame wavelength and the fact that it is likely produced in 
the inner, probably most obscured part of the NLR, it may suffer 
from extinction. While values above few hundreds are definitively 
an indication of lack of obscuration, almost all objects with ratios 
below 15 are Compton thick (see the luminosity plot in Fig. 2 of 
Gilli et al. 2010 for objects at comparable redshift). Hence, no firm 
conclusion on the amount of obscuration close to 109104 can be 
inferred from the X-ray/[Ne V] ratio. 



6 SUMMARY 

We have presented sensitive Chandra data (i.e., with significantly 
higher photon statistics than those reported by 101) and UV/optical 
to far-IR SED fitting of the narrow-line quasar 109104 at z = 0.442 
in order to characterize the properties of the matter responsible for 
the X-ray absorption and reprocessing of the nuclear radiation in 
the mid-IR (i.e., the torus envisaged by AGN unification schemes). 
The source detection in the PDS instrument onboard BeppoSAX 
more than ten years ago has brought 109104 to be classified as a 
Compton-thick AGN (F00). 

In the following, we summarize the main results achieved with 
the multi-wavelength analysis presented in this work. 

• The X-ray spectrum is in favour of a Compton-thin hypothesis, 
which has also been suggested by P07 from the analysis of XMM- 
Newton data. 

• The 54-month Swift BAT map in the 15-30 keV energy range 
has clearly shown that excess hard X-ray emission, if present (since 
formally below the BAT detection threshold adopted for source cat- 
alogs), is likely to be ascribed to NGC 2785, a nearby (z = 0.009) 
Type 2 AGN at ss 17' from 109104. The properties of its SDSS 
spectrum are suggestive of heavy extinction. 

• Further support to the Compton-thin hypothesis comes from 
observations at longer wavelengths, where the obscuring matter 
(torus), responsible for all of the observed emission up to 60iim, 
can be well parameterized, providing a x(9.7)=3.0, consistent with 
the absorption measured through X-rays (using a Galactic gas- 
to-dust ratio conversion). The covering factor of the absorber is 
fa 90 per cent. Broad-band photometric data are well reproduced by 
a stellar component plus a smooth dust distribution (as described 
by F06); the silicate absorption feature at 9.7(im observed in the 
Spitzer IRS spectrum is also properly accounted for by the same 
torus model. On the other hand, clumpy-torus models (N08a,b), 
which likely provide a more realistic (physical) description of the 
obscuring matter (e.g., Tristram et al. 2007), show problems in 
reproducing the source photometry and spectroscopy at the same 
time. 

• Using the mid-IR (at 12.3fi.rn and 5.8(im, derived directly from 
the IRS spectrum) and the [O III] emission to predict the hard X-ray 
luminosity of 109104 provides some space for the presence of addi- 
tional absorption to be added to the value of column density derived 
from Chandra data. 

The scientific case of 109104 strongly supports the importance 
of broad-band, possibly simultaneous X-ray spectral coverage to 
define the source complexities properly and provide a robust esti- 
mate of the amount of absorption towards the source. In this regard, 
hard X-ray spectroscopic and imaging capabilities, as those pro- 



vided by NuSTAR (planned to be launched in 2012), might provide 
the final answer on the nature of the absorber in 109104. 
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